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HYDROGRAPH ANALYSIS BY THE METHOD OF LEAST SQUARES 


Willard M. Snyder,! A.M. ASCE 


ABSTRACT 


The hydrographs of ten storms were analyzed by a technique based on the 
method of least squares. By iterative solution coefficients of distribution of 
total runoff were obtained simultaneously with estimates of the total runoff 
volume. 

Five of the storms were taken from a watershed with changing cover con- 
dition. The other five were taken from a smaller watershed with stable cover 
condition. The coefficients of runoff distribution for the storms on changing 
cover show a change in shape of the total storm hydrograph. The distribution 
coefficients for the storms on stable cover do not show a change. 


INTRODUCTION 


The analysis of the hydrology of small watersheds is usually based, in part, 
upon storm lists which are used to relate volume of run-off and peak discharge 
to various parameters of rainfall and antecedent conditions. In addition, in- 
dividual storms are selected from the storm lists to obtain characteristic 
hydrographs of unit volumes of runoff, and to obtain average rates of reces- 
sion from peak flows. 

Many valuable conclusions regarding the effect of land management upon 
hydrologic ar be drawn from these studies. However, there 
are certain limitations.\!,2) surface runoff, for example, conceivably could 
be delayed by improved cover. However, if no change took place in the vol- 
ume of surface or ground-water runoff, volumetric studies will not detect 
changes in the hydrograph. A study of peak discharges could show that a 
change had been caused by land management. However, this change in peak 
discharge is only one aspect of the total change in shape of the hydrograph. 

Studies of runoff and peak discharge may be supplemented by a study of 
the hydrographs of selected storms. Here another difficulty is sometimes 
encountered. To compare storm hydrographs the storms must occur under 
similar antecedent conditions and must have similar rainfall amounts and in- 
tensity patterns. Sets of storms meeting these requirements are difficult to 
find even when the hydrologic record is long. 

The largest storms in a hydrologic record may be so complex that they 
cannot readily be used in the hydrologic analysis. To compute the volume of 
surface runoff of these complex storms serious assumptions must be made 
concerning the separation of surface and ground-water runoff components. 
Such storms can hardly ever be used in direct hydrograph comparisons be- 
cause the complex rainfall pattern has little chance of being duplicated in 
another storm. Yet, it is these largest storms of record which are often the 
most important from a hydrologic viewpoint. 


1. Hydr. Engr. Tennessee Valley Authority, Hydraulic Data Branch, Hydr. 
Section, Knoxville, Tenn. 


793-1 


A technique by which complex storm hydrographs can be analyzed by the 
method of least squares is presented in this paper. The basic technique is not 
new.(3) However, it has been extended to the development of approximate unit 
hydrographs of total storm runoff. The technique applies to storms with sim- 
ple rainfall patterns, as well as to storms with complex rainfall patterns. 

Only winter-time storms of fairly long duration were used in the development. 


The Watersheds 


This technique of hydrograph analysis by the method of least squares was 
developed using rainfall and runoff data observed by the Hydraulic Data Branch 
of the Tennessee Valley Authority in the White Hollow Experimental Watershed 
in Union County, Tennessee. The method was tested and refined using data 
from Watershed 3 of the Tennessee Valley Authority-North Carolina State 
College Cooperative Research Project at Waynesville, North Carolina. 

The White Hollow Watershed has been described in detail in an earlier pub- 
lication. 4) Briefly, the watershed is an area of 1,715 acres which was ac- 
quired by the Tennessee Valley Authority in connection with its land acquisi- 
tion program for Norris Reservoir. Some 30 families who resided there in 
1934 were moved out of the watershed, and the area was taken out of cultiva- 
tion. In 1934, about two-thirds of the watershed was forest which had been 
subject to burning, grazing, and cutting, and the remainder was open. 

When the watershed was taken out of cultivation it was placed under pro- 
tection. Extensive erosion control work was carried out during 1934 and 1935. 
Tree planting in open and forested areas was carried on in 1934 and 1935, and 
during the period 1938 and 1942. A dense weed growth became established on 
the formerly cultivated areas and these gradually gave way to herbaceous and 
woody species. By 1946, virtually all of the area was classified as forest. 

The few remaining open areas were homesites of less than one acre. With 
continued protection and management, the protective cover of the watershed 
has continued to strengthen. 

Watershed 3 of the Cooperative Research Project in Western North 
Carolina is an area of 3.5 acres, located about 12 miles north of Waynesville, 
North Carolina. The watershed has an elliptical shape, being about 550 feet 
long and 350 feet wide. The slopes are rather uniform and moderately steep, 
averaging about 26 percent. 

The soil of Watershed 3 is deep to bedrock. The area, once badly eroded, 
had been improved by the owner over a period of years as a test demonstra- 
tion farmer. The ground surface is covered by an improved pasture. 

Runoff from the watershed is measured by a 2-foot type H flume equipped 
with a recording gage. A cut-off wall of sheet metal extends from the flume 
to a rock ledge about six to eight feet below the surface, so that virtually all 
runoff from Watershed 3 is measured. A recording rain gage is located in 
the area. 

There is no continuous base flow from Watershed 3. Runoff occurs for 
comparatively brief periods after intense or prolonged rain. 


The Method of Analysis 


When a unit hydrograph of a defined duration is available for a watershed, 
it can be used to synthesize a complex storm hydrograph resulting from suc- 
cessive periods of runoff-producing rainfall. To do this, the storm is divided 
into periods equal to the defined duration of the unit hydrograph. The volume 
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of runoff in each of these periods is distributed by multiplying by the appro- 
priate coefficients of the unit hydrograph. Each hydrograph ordinate is then 
the sum of the separate products of runoff times unit hydrograph coefficient. 
Computation of the ordinates of a storm hydrograph may be expressed by 
a linear equation as follows: 
Ye = U,RO, U,RO, + U,RO, (1) 
where Ys is the hydrograph ordinate at the end of 5 periods of runoff, 
RO, to ROs are the amounts of runoff during the successive periods, and 
Ua, to Up are the first 5 coefficients or ordinates of the unit hydrograph. 


The ordinate of the storm hydrograph after 5 periods of runoff was used as 
an example. A similar equation can be written for any ordinate of the hydro- 
graph. 

By convention, several assumptions and interpretations are implied in 
Equation 1. The coefficients U are usually distribution parameters of surface 
runoff, and all the coefficients U cover equal periods of time. The RO values 
represent surface runoff volumes during successive equal periods of times. 
Short periods and a large number of coefficients U may be required to define 
adequately the entire runoff distribution from a small watershed having a 
sharp-crested hydrograph. 

For the analysis presented in this paper the form of Equation 1 was retained 
but some definitions were changed. The coefficients U were assigned effective 
distribution times as shown in Table 1. The total storm hydrograph was used, 
instead of the hydrograph of surface runoff. That is, the coefficients U are so 
defined that a distribution of total runoff, rather than surface runoff, is accom- 
plished. The representation of the hydrograph of total runoff by the distribu- 
tion coefficients U is shown in Figure 1. In this figure a continuous hydro- 
graph has been sketched through a plotted set of coefficients. The coefficients 
are average ordinates of the hydrograph for the time intervals defined in 
Table 1. 

Using the distribution coefficients for the White Hollow Watershed, the 
equations for successive hourly ordinates of a storm hydrograph are: 


U,RO, 
U, (RO, + RO, ) 
= U, (R05 + RO, ) + 


U, (RO), + R03) + U,(R0, + RO,) 


U, (RO, + RO) ) + Up (RO, + RO, ) + 


U, (RO, + ROg) + Up (RO), + + U, (RO, 


= J, 
U + + UpRO, 


The storm ordinates Y in Equations 2 are defined as the increase in flow 
during the storm over the flow antecedent to the storm. As such the increase 
includes all surface, subsurface, and ground-water flow associated with the 
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current storm. Equations 2 can be extended to cover a total hydrograph 
length of 64 hours from the beginning of runoff. The runoff producing rainfall 
may be any duration up to 64 hours. Durations longer than 64 hours can be 
analyzed if the runoff is considered to have no further effect on the hydrograph 
64 hours after the runoff occurred. This assumption will introduce some er- 
ror, mainly on the recession side of the hydrograph. 

It was stated in the discussion preceding Equations 2 that each ordinate of 
a storm hydrograph is the sum of the products of runoff volumes times dis- 
tribution coefficients. In Equations 2, it can be seen that the product terms 
containing Ua, are the portions of the hydrograph ordinates caused by runoff 
occurring from 0 to 2 hours prior to the ordinates. The terms containing Up 
are the portions of the ordinates caused by runoff occurring from 2 to 4 hours 
prior to the ordinates. The terms containing U,, therefore, measure that 
portion of the watershed total runoff which has a time-of-travel to the gaging 
station varying from 0 to 2 hours, The terms containing Up, measure that 
portion of the watershed total runoff which has a time-of-travel from 2 to 4 
hours. Each coefficient may be regarded as measuring the portion of total 
runoff with a time-of-travel corresponding to the effective times given in 
Table 1. 

It is believed that description of the total storm runoff by time-of-travel is 
a realistic approach. For example, surface runoff is not often observed in the 
White Hollow Watershed except from the few access roads and along the road 
ditches. Where these ditches enter the watershed proper this surface water 
usually infiltrates and does not reach the stream system as overland flow. 
Conversely, it has frequently been observed that near the valley bottoms there 
is effluent water from subsurface flow which travels a short distance as over- 
land flow to the stream system. There is apparently no very real separation 
of runoff into fixed components of surface, subsurface, and ground-water 
runoff. 

The technique of analysis to be presented requires the application of 
Equations 2 to the total storm hydrographs of selected storms by the method 
of least squares. A set of simultaneous equations must be solved, and the 
number of equations in the set equals the number of coefficients U. Sets of 
more than about six to eight simultaneous equations are difficult and laborious 
of solution with a desk calculator. 

In actual practice, runoff by hourly amounts, as indicated in Equations 2 is 
not accurately known. The normal practice in bee a hydrograph of sur- 
face runoff is to estimate these hourly runoff amounts.(°) In this analysis, 
however, a method is presented for obtaining the hourly runoff amounts by 
iteration. The coefficients U and the runoff amounts RO are obtained by suc- 
cessive approximations. 

For the first approximation hourly rainfall amounts instead of hourly run- 
off amounts were used to evaluate a set of coefficients. Equations 2 were 
revised as shown in Equations 3. 


UR, 


U.(R, +R 
+ * Ry) 


+ (8, + Ro) + 


at 
(3) 
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In Equations 3, the successive ordinates of the storm hydrograph were 
predicted by successive products of rainfall, R, and distribution coefficients, 
U. In addition, a constant term, a, was added to each equation to allow some 
compensation for the difference between rainfall and runoff. 

The first step in analysis consisted of fitting Equations 3 to observed data 
by multiple regression. For a specific storm, the observed values of the de- 
pendent variable, y, are the successive hourly ordinates of the storm hydro- 
graph. The observed values of the independent variables are the hourly rain- 
fall amounts, distributed and combined as indicated by the parenthesis in 
Equations 3. The coefficients, U, and the constant, a, are the regression co- 
efficients estimated by the least squares solution. 

Table 2 shows the rainfall and discharge data for the storm of January 1, 
1937, on the White Hollow Watershed. These data have been tabulated in the 
form required for multiple regression. Rainfall is tabulated by hourly 
amounts. These amounts are combined by periods corresponding to the time- 
definition of the coefficients. The combined amounts are tabulated under col- 
umn headings Ra, Rp, etc., to show the association with the respective co- 
efficients. The ordinates of the hydrograph, y, are also tabulated. These or- 
dinates were obtained by projecting the recession curve of the flow antecedent 
to the storm. This recession curve was subtracted from the total observed 
hydrograph to obtain the hydrograph of the current storm. 

A discussion of the application of the method of least squares to the date in 
Table 2 will not be given here. Detailed descriptions of this method may be 
found in any standard Statistical text. 

The coefficients U resulting from the analysis of data in Table 2 are shown 
in Figure 1. The coefficients are shown in proper time relationship to one 
hour of runoff-producing rainfall. It is evident from Figure 1 that the coeffi- 
cients of Equations 3 define an approximate unit hydrograph. It should be 
pointed out that the coefficients in Figure 1 have not been adjusted. Therefore, 
the area bounded by the coefficients does not necessarily represent one inch 
of total runoff from the watershed. 

The coefficients shown in Figure 1 can be used to predict the storm hydro- 
graph from the rainfall data in Table 1. This procedure is almost identical 
to the usual method of predicting a hydrograph when the runoff amounts and 
the unit hydrograph are known. The constant, a, in Equations 3 also enters 
as an additive term to each ordinate prediction. The hydrograph predicted by 
this first solution and the observed hydrograph for the storm of January 1, 
1937, are shown in Figure 2. 

There is generally fair agreement between the observed and predicted 
hydrographs. However, the predicted hydrograph is too high for about the 
first six hours of rise, and also falls too rapidly during the recession portion 
of the hydrograph. These errors of fit have a rational explanation. The basic 
failure of Equations 3 lies in the inadequacy of the constant, a, to express the 
changes in the difference between rainfall and runoff during the storm. This 
difference, or loss rate, is high at the beginning of a storm. During the 
course of the storm, as interception and retention losses are satisfied, and as 
infiltration capacity decreases, loss rates become smaller and approach a 
constant value for long-duration storms. 

In the storm of January 1, 1937, the first burst of rainfall produced very 
little runoff. As rainfall intensities again increased, rates of runoff increased 
greatly. The constant, a, in Equations 3 obviously cannot express this chang- 
ing relationship of rainfall to runoff. 

The failure of the hydrograph predicting equation during the recession is 
a consequence of the method of least squares. The sum of the residual errors 
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of prediction must be zero. The over-prediction at the beginning of the storm 
must be balanced by under-prediction during the recession. 

If allowance is made for the inability of Equations 3 to express the chang- 
ing loss rate, the general fit of the predicted storm hydrograph to the observed 
storm hydrograph is good. The coefficients U, therefore, probably are a rea- 
sonable first estimate of a unit hydrograph. 

When the unit hydrograph is known, the loss for each period of rainfall dur- 
ing a storm can be estimated by the following method. Loss may be expressed 
by 


(4) 


As expressed in Equation 4, L, the loss is defined as all rainfall which does 
not become runoff. This loss must be negative or zero. It cannot be positive, 
since runoff cannot be greater than rainfall. 

Each predicted ordinate of the hydrograph, as computed by Equations 3, has 
associated with it some error of prediction, E. This error may be defined as 
the observed ordinate minus the predicted ordinate, or, in equation form 


a - - 


In a) 


where y is the observed ordinate and (Y' - a) is the predicted ordinate, un- 
adjusted for the constant a. 

In Equation 2 the computed ordinate, Y, is an estimate of the observed or- 
dinate, y. The various terms, URO, from Equation 2 can then be substituted 
for y. Similarly, from Equation 3, it can be seen that the terms, UR, can be 
substituted for Y' - a. Carrying out both these substitutions produces 


B, U,RO, UR, 
U, (RO, + RO, ) - U, + 


2 
U, (Ro, + + - + Ro) - 


By rearranging terms and substituting L = RO - R from Equation 4 there 
results 
E, * U,(L, + L,) 


(5) 


(4) 
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In Equation 5-1, estimated values of E; and U, are known from the mul- 
tiple regression analysis, and the equation can be solved for Lj. This value 
of Ly and the estimated value of Ep allows solution of Equation 5-2 for Lo. 
Values of L can thus be computed for each predicted ordinate Y'. Obviously, 
values of L = 0 should be used for those periods during which no rain falls. 

The values for L will be negative near the beginning of a storm and positive 
near the end of the storm, since the hydrograph was over-predicted at the 
beginning, and under-predicted toward the end. Physically, the values for L 
can only vary between the limits of zero and - R, and these limits should be 
retained during computation. In effect, values of L were computed for all 
periods during which rain occurred. Only negative values of L were retained, 
and these values were not allowed to be larger (negative) than the rainfall for 
the period. Positive values of L were set equal to zero. 

The reliability of these first estimates of loss depends upon the first esti- 
mates of the coefficients, U, and upon the computed errors, E. The estimate 
of E is probably less reliable than the estimate of U, since E contains the 
statistically random fluctuations in the rainfall-hydrograph correlation. 
Therefore, at this point in the analysis the estimated losses should only be 
regarded as an analytical device. 

Following the procedure outlined above, the losses were computed for the 
storm of January 1, 1937, on the White Hollow Watershed. These losses are 
shown in Table 3, together with the solution of Equation 4 to obtain hourly val- 
ues of runoff, RO. The right-hand side of Table 3 shows the values of runoff 
combined for association with the respective distribution coefficients. 

The next step in the analysis was to fit Equations 2 to the data in Table 3 
by the method of least squares. This fitting is used to obtain a new estimate 
of the distribution coefficients, U. This second set of coefficients is shown 
with the first set in Table 4. 

The second set of coefficients was used in Equations 2 to again predict the 
storm hydrograph. This second prediction is shown with the first prediction 
and the observed hydrograph in Figure 2. The hourly runoff values from 
Table 3 have also been plotted in Figure 2. The agreement between the sec- 
ond prediction and the observed hydrograph is very good. Some additional 
adjustments could be made near the crest of the hydrograph, and also follow- 
ing the small burst of rain which occurred during the recession. 

Since there are only small differences in the two sets of coefficients in 
Table 4, and since the predicted storm hydrograph agrees well with the ob- 
served storm hydrograph, the analysis for the storm of January 1, 1937, was 
stopped at this point. However, it should be pointed out that the second set of 
distribution coefficients could have been used to re-evaluate the estimated 
losses, a new set of hourly runoff values could have been obtained, and from 
these a third estimate of the distribution coefficients could have been made. 
This iteration was performed on several of the other storms presented in this 
paper. The number of trials needed depends upon the desired precision of 
the set of distribution coefficients. 


Refinement of the Method of Analysis 


The method of analysis presented above was used to analyze five storms 
on the White Hollow Watershed. The method was later improved by a modifi- 
cation of the method of computing the losses. The revised method was used 
to analyze five storms on Waynesville Watershed 3. 
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In the initial determination of a set of distribution coefficients by Equations 
3, the constant, a, was included to allow some expression of loss. The initial 
solution can be improved by using a “dummy” variable to approximate a chang- 
ing loss rate. This “dummy” variable replaced the constant, a, in Equations 
3, by aterm apS,. If 8, be assigned values of unity for specific ordinates of 
the hydrograph, the solution by least squares will produce values for a,, a 
regression coefficient, which apply only to the specified ordinates. For ex- 
ample, Equations 3 might be re-written as: 


J). 
Y,' 4,5) 
Y,' a8, + + R,) 

R 4 ULF 


The solution of Equations 3’ is controlled by setting 8, equal to unity for 
ordinate Y;', and setting it equal to zero for Yo’ and Y3'. Similarly, So is 
set equal to unity for ordinate Y»' , and set equal to zero for ordinates Y;' 
and Y;'. By this “on-off” variation of the dummy variable, S,, the coefficients 
ap, are associated only with the specified ordinates. This coefficient, a,, 
which is an approximation of a changing loss-rate, can be used in place of the 
constant, a, in the development of Equations 5. 

In practice about four values of a, were determined: a, applying to the 
first quarter of the storm, ag to the second quarter, and soon. Because of 
the large number of zeros in the simultaneous normal equations when an “on- 
off” technique is used, the number of equations can be increased to about 10 
or 12, and can still be solved using a desk calculator. 

In addition to the use of a dummy variable, the computation of losses was 
further modified by a change in the solution of Equations 5. Step-wise solution 
of these equations is unsatisfactory if poor initial estimates of E, and U, 
produce a poor estimate of L,, and consequently poor estimates of all the fol- 
lowing values of L. By a technique based on least squares the errors in the 
values of L may be minimized simultaneously over several rainfall periods. 

Equations 5 may be expanded and rearranged as follows: 


E, = L,U, + LU, 


= T 
LU, + + 


= LU, + Lt, + 1,0, + + Lely (5a) 
Be + + LU, + LU, LeU, + 
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When the terms are arranged as shown in Equations 5a it can be seen that 
the values of L can be computed from estimates of the coefficients U and er- 
rors E by multiple regression. If the losses for the first six hours of runoff- 
producing rainfall are computed simultaneously by the method of least squares, 
a better estimate of these losses is obtained than if they are computed by pro- 
gressive single solution of Equation 5. 

Due to the limitations of a desk calculator, only about six to eight values of 
L can be computed simultaneously. It is necessary to use these values to com- 
pute the next group of six to eight values of L. In effect, the values of L are 
computed in sets of six to eight, for as many sets as necessary to cover the 
entire storm period to that point where the values of L become positive. 

The method of successive calculations is indicated below. The second set 
of loss equations is: 


TT T 1 
= LU, + LU, + LeU, + 


+ 


+ 


LeU, + LeU, 


+ 


+ + Lee + 
+ LU Lgl, + Lou, + + 
The prediction errors Eq through Ej9 can be adjusted for the effect of the 
first set of six losses by setting 


Eg Lu, L, Uo LU, Le LU, 


Substituting these values in Equations 5b there results: 


1's ] 


U, U 

ra” Lg A 

Lgl, Lgl, LU, 


Bo ' 
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+ Lgl, 

Ey + + 13% + LU + 
+ Lgl, + LU, 

* + LU, + + LU, + (5b) 


Note that the right-hand sides of Equations 5a and 5c are identical in the 
arrangement of the coefficients U. This greatly reduces the amount of compu- 
tation required in estimating the second set of losses. In Equations 5a the co- 
efficients U are considered the independent variables in the solution by least 
squares. Since the coefficients U are identical in Equations 5a and 5c the si- 
multaneous normal equations for solution of these equations are identical in 
those portions that obtain from the independent variables. Therefore, only 
that portion of the solution by least squares which is altered by introduction 

of a new dependent variable need be re-solved for the second set of losses. 

This process of analysis can, of course, be carried forward to the solution 
of a third set of losses. As the losses are computed, it should be remembered 
that they cannot be greater (negative) than the rainfall, nor less than zero. 

Any computed loss values falling outside these limits should be set equal to 
the rainfall, or to zero, as appropriate, before the losses are used in comput- 
ing the next set. 

After all loss values have been estimated, runoff values can be estimated 
by subtracting the loss from rainfall. From this point the hydrograph analysis 
follows the iterative procedure described earlier, except that any re-evalua- 
tion of losses is accomplished by solution of Equations 5a and 5c by the meth- 
od of least squares. The desired precision again determines the number of 
re-correlations. 


Results of Analysis 


White Hollow Watershed 


Four storms were analyzed in addition to the storm of January 1, 1937, 
which was used as an illustration. Data for the five storms are shown in 
Table 5. The rainfall amounts are the total storm amounts. The storm losses 
are the summations of the L values, and the runoff amounts are the summa- 
tions of the RO values, as shown in Table 3. 

The rainfall hyetographs, the stream flow hydrographs, and the results of 
analysis for the storms of January 1937 and February 1945 are shown in 
Figures 3 ani 4. Time-of-travel separations of the hydrographs are also 
shown in these figures. The lines of separation were computei by dropping 
terms from the right-hand side of Equations 2 after computing the ordinates 
of the predicted hydrograph. Thus, the line labeled “2 hours” is the summa- 
tion of all terms in Equations 2 except U,. The line labeled “4 hours” is the 
summation of all terms except U, and ese lines of separation in Fig- 
ures 3 and 4 are very similar in appearance/to conventional lines of flow 
separation in hydrograph analysis. It sh be emphasized, however, that 
these lines result from the method of analygis. They are not arbitrary lines 
of separation which must be used as a starting point in the analysis of a hydro- 
graph of surface runoff. ) 

The final sets of coefficients U for the/five storms are shown in Table 6. 
The coefficients as shown here have been adjusted by simple proportion so 
that the total flow represented by the coefficients equals one inch of runoff 
from the watershed. The factors required for this proportional adjustment 
are included in Table 6. 
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The storm of January 30, 1950, was of long duration. Consequently, the 
definitions of effective times for the distribution coefficients was changed for 
this storm. The total time covered by the coefficients was 128 hours. Only 
that portion of the set of coefficients equal to the first 64 hours is shown in 
Table 6. The factor of 1.15 to convert to unit runoff applies to the 64-hour 
portion of the distribution coefficients, so that these adjusted coefficients are 
comparable to those for the other four storms. 

The distribution coefficients for total runoff listed in Table 6 show a con- 
sistent shift. The storms of January 1937 and February 1939 had highest val- 
ues for the coefficients U, and U,; that is, had peak distribution of runoff 
from 0 to 4 hours after rdinfall. e three later storms show that the peak 
distribution of runoff occurred about 4 to 8 hours after rainfall. The coeffi- 
cients for those runoff times longer than about 8 hours increased from the 
earlier to the later storms, showing that runoff time-of-travel increased. 

It might be thought that the change in runoff distribution is only apparent, 
caused by the progressively smaller storm rainfall amounts for the first four 
storms. Yet the storm with the largest rainfall amount, 5.05 inches, was the 
latest storm, January, 1950. This storm has distribution coefficients very 
similar to those for the smaller storms of February 1945, and March 1949. 

The differences in distribution coefficients does not appear to be primarily 
the result of differences in rainfall intensity. The storm of February 1939 
had the highest rainfall intensity, and this intense rain occurred after some 
immediately prior rain which would have tended to lower the infiltration rate. 
Yet this storm shows some delay in runoff distribution when the coefficients 
are compared with those for the earlier storm of January 1937. 


Waynesville Watershed 3 


The hydrographs of five storms were analyzed. Data for the five storms 
are shown in Table 7. Only that rainfall which occurred during the period of 
runoff was used in the hydrograph-analysis. Table 7 shows those rainfall 
amounts which were “lost” before runoff began. Some error is involved in 
this procedure, but it is believed the error is small. It is necessary to make 
some adjustment of this kind on Watershed 3 where 10-minute rainfall incre- 
ments were used. A large number of these short periods might be associated 
with zero ordinates before the beginning of runoff and produce an unrealistic 
correlation. Also, since there is usually no base flow from Watershed 3, 
initial losses of detention and retention tend to be high. On the larger White 
Hollow Watershed, where hourly rainfall amounts were used, usually only one 
or two periods of rainfall occur before some runoff begins, and the handling of 
data is then not critical. 

The hydrographs for the January 1950 and February 1953 storms are 
plotted in Figures 5 and 6. Rainfall is shown by 10-minute increments. The 
predicted hydrographs and the 10-minute runoff increments are included. The 
runoff distribution coefficients resulting from the hydrograph analyses are 
shown in Table 8. These coefficients have been adjusted to equal one inch of 
total runoff from the watershed. 

The distribution coefficients shown in Table 8 are fairly stable in arrange- 
ment. The coefficients U,, to U_, particularly, show little variation, thus 
demonstrating that the diftributfon of runoff beyond about 2-1/2 hours time- 
of-travel has little variation. The peak of the runoff distribution is either at 
U, or Uc, or, for the storm of January 1954, a double peak at U, and Uc. 
The five storms which were analyzed were all different in their basic rainfall 
intensity patterns. Yet only two patterns of distribution of runoff are evident 
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in Table 8, and there is no general shift of the distribution coefficients with 
time. 

A close examination of the rainfall patterns and hydrographs for the five 
storms showed there were instances where high rainfall intensities seemed to 
produce almost immediate response in the hydrograph. It may be possible 
that the peak coefficient U, is associated with this rapid response of the hydro- 


graph. 


Comparison of Results 


The change in the distribution of runoff from White Hollow Watershed be- 
comes more evident when the distribution coefficients in Tables 6 and 8 are 
used to compare the time required for equivalent proportions of runoff. The 
times for 25 percent and 50 percent of total runoff were computed for the ten 
storms. These times are shown plotted agains the year of storm occurrence 
in Figure 7. The time required for 25 percent of the total runoff from White 
Hollow Watershed increased from about 4 hours in the 1937 and 1939 storms 
to more than 8 hours in the 1945, 1949, and 1950 storms. The time for 50 per- 
cent of total runoff increased from 9 hours in the two earliest storms to more 
than 15 hours in the three later storms. 

There is no time-trend in the time required for 25 percent or 50 percent of 
the total runoff from Waynesville Watershed 3. The major variations in the 
time-distribution of runoff were caused by the single peak distribution of run- 
off at coefficient Uc for the storms in 1951 and 1952. 


CONCLUSIONS 


The hydrographs of total runoff of complex storms can be analyzed by the 
method of least squares. Sets of coefficients which define the time-distribu- 
tion of runoff result from the analysis. 

By an iterative procedure the rainfall loss, and, therefore, an estimate of 
volume of runoff, can be determined simultaneously with the distribution co- 
efficients. 

The distribution coefficients for five storms on White Hollow Watershed, 
an area of 1715 acres, show that winter-time runoff was delayed as the forest 
cover improved. 

The coefficients for five storms on Waynesville Watershed 3, an area of 
3.5 acres, show no major change in the time distribution of runoff. This 
watershed was under a stable cover of improved pasture. 
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TABLE 1 


TIME DEFINITION OF DISTRIBUTION COEFFICIENTS 


Waynesville 
Coefficient White Hollow Watershed 3 


Effective time from Effective time from 
beginning of 1 hour beginning of 10 minutes 
of runoff in hours of runoff in 10-minute 
Periods 
Oto 2 Oto 
2to ite 3 
k to 8 3 to 7 
8 to 16 7 to 15 
16 to 32 15 to 31 
32 to Gy 31 to 63 


Not Used 63 to127 


Uy 
UB 
Up 
Up 
Up 
Ug 
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TABLE 2 
HYDROGRAPH ANALYSIS 
WHITE HOLLOW WATERSHED 


Storm of January 1, 1937 
Data for Solution of Equations 3 
All Data in Inches Except as Noted 


ax 
a 


wr 
ar 


10 
2k 
15 
12 
03 
06 
207 
1 
5 
3 
7 
8 
6 
8 
7 
3 
1 
1 


0.1 
0.3 
0.2 
0.3 
0.1 
0.1 
0.3 
0.3 
0.1 
0.1 
0.0 
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Rain- 
Date Time fall Ry R, Ry Rp | 
1 6p -10 
7p 
8p 039 10 
9p 027 
10p 15 .39 
llp .27 
Mdt 
2 la | 
Sa oan 61 89 
$ 6a 39 99 
7a 111 
8a 83 e777 129 
10a 50 .54 1.12 .93 10 132 
lp 0 OL .2k 1.09 1.3% .61 91 
2p 0 0 1.0h 1.51 .64 80 
3p 0 0 -O1 -99 1.61 -70 bz! 
lp 0.08 -08 O -62 1.92 off 64 
Sp 0.06 O -25 2.18 59 
6p 0.02 .08 .08 .12 2.16 1.03 55 
7p 0.08 1.93 1.37 52 
8p 0.06 .08 -08 1.71 1.60 
9p 0 06 1.34 1.97 L6 
10p 0 lh .16 1.16 2.15 L6 
llp .24 1.00 2.31 L6 
Mdt 0 222 -70 2.69 Lu 
3 6a 30 2.67 Gh Fi 
Noon 2. 1.60 21.2 
6p 3.19 16.3 
Mdt 3.39 12.4 
4 6a 0 3.61 9.5 
Noon 3.12 7.6 
Total 3.61 


TABLE 3 
HYDROGRAPH ANALYSIS 
WHITE HOLLOW WATERSHED 


Storm of January 1, 1937 
Data for Solution of Equations 2 
All Data in Inches Except as Noted 


RO RO RO, RO 


A RO, RO, 


Total 


@ 
WNW 


0 
0 
0 
0 
0 
0 
0 
0 
0. 
0 
0 
0 
0 
0 
0 
6) 
0 


OA 


. 
wo 


ERS 
Or 


0900 0 


« 


oo 


Re 
oO 


3.61 1.111 2.499 


-002 
-016 
-028 
+037 
079 
2129 
2239 
+769 
1.159 
1.339 
1.499 
1.879 
2.L62 
1.710 
0.12 
0 


-037 
-789 
2.379 
2.499 
2.499 
2.483 


a 


On 


+ 
1 6p 0.002 .002 
7p 0.005 .907 
8p 0.009 .002 
op 0.012 .021 .007 
10p 0.009 .021 .002 
llp 0.017 .026 .021 .007 
Mdt 0.025 .O42 .021 .016 
2 la -060 9.050 .075 .026 .028 
2a 0.110 .160 .Ol2 .035 .002 
3a 0.320 .430 .075 .047 .007 
ha 0.23 .55 .160 .063 .016 76 
Sa 0.37 .60 .430 .101 .028 89 
6a 0.18 .55 .55 .202 .037 99 
Ta 0.16 .34 .60 .505 .054 111 
Ba 0.38 .55 .710 .079 129 
9a «3h: 1.08. 133 
10a «50 «Sh 132 
lla Onl] .75 .552 116 
Noon 0.01 .12 .501.09 .773 103 
lp oO 1.09 1.131 91 
2p 0 0 -12 1.04 1.302 80 
3p -01 0.99 1.4L5 71 
lp 0.08 0. 0 -62 1.600 64 
Sp 0.06 0. 59 
6p 0.02 0. oie 55 
7p 0.08 O. -O1 1.93 52 
8p 0.06 0. 1.71 L9 
9p O 0 1.3L L6 
10p 1.16 
llp 1.00 Li6 
Mdt 0.62 Lh 
3 6a 0 31 
Noon 21.2 
6p 16.3 
Mdt 12.4 
L 6a 9.5 
Noon 7.6 
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TABLE 
DISTRIBUTION COEFFICIENTS FOR TWO TRIALS 


WHITE HOLLOW WATERSHED 


Storm of January 1, 1937 


First 
Coefficient Trial 


92 


B 63 


U 


C 


TABLE 5 
STORM DATA 
WHITE HOLLOW WATERSHED 


All Data in Inches 


Storm Date Kainfall Runoff 
January 1, 1937 3.61 2.50 
February 2, 1939 2.79 
February 17, 1945 
March 17, 1949 1.87 


January 29, 1950 5.05 


Trial 
76 
37 
Up 26 23 
Ug 17 13 
Up 6 4 
1.54 3.51 
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TABLE 6 
RUNOFF DISTRIBUTION COEFFICIENTS 
WHITE HOLLOW WATERSHED 
Effective Jan. Feb. Feb. Effective Jan. 


Coefficient Time 1937 1939 1945 Time 1950 
Hours Hours 


O- 2 0 - 
2-4 3 
6 
6 - 16 10 


TABLE 7 
STORM DATA 
WAYNESVILLE WATERSHED 3 
All Data in Inches 


Rainfall Loss 
Total Initial During Period During Period 7 
Storm Date Rainfall Loss of Runoff Runoff of Runoff 
January 19, 1950 1.29 
March 6, 1951 1.06 


January 9, 1.92 


February 20, 195 2.00 0.469 


January 20, 1! 2.06 0.64 


Up 6 64 
Ue 10 66 
\ 14 62 
Up 16 = 32 23 18 29 26 14 - 22 ko 
32 = 64 7 8 16 22 17 
ho 64 13 
Factor 1.69 0.66 1.33 1.62 1.15 
otal 
Loss 
.68 
62 
(0.26 1.46 1.19 0.47 0.73 
ee 1.31 0.90 0.41 1.10 
ee 1.42 0.56 0.56 1.20 
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TABLE 8 
RUNOFF DISTRIBUTION COEFFICIENTS 
WAYNESVILLE WATERSHED 3 


Effective Jan. Mar. 
Coefficient Time 1950 1951 
Minutes 


O- 10 3-55 0.17 


10 - 30 0.49 


Factor 


Jan. Feb. Jan. 
1952 1953 1954 
Up 1.23 1.06 
Uy 30 = 70 6.96 
Ug 150 - 310 0,19 0.34 0.109 G23. Ge 
U; 630 -1270 0.0; 0.066 0.0 0.02 0 
149 2.00 3.96 2:0h 2.25 
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